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Abstract 

■ 

In this paper we discuss the representation of the joint probability density function of perfectly correlated continuous random 
variables, i.e., with correlation coefficients p=±l, by Dirac's 5-function. We also show how this representation allows to define 
Dirac's <5-function as the ratio between bivariate distributions and the marginal distribution in the limit p — > ±1, whenever this 
limit exists. We illustrate this with the example of the bivariate Rice distribution. 

I. Introduction 

^ \ The performance evaluation of wireless communications systems relies on the analysis of the behavior of signals modeled 
as random variables or random processes that obey prescribed probability disnibution laws. The Pearson product-moment 
correlation coefficient is commonly used to quantify the strength of the relationship between two random variables. For 
example, in wireless communications, two-antenna diversity systems take advantage of low signal correlation. The diversity 
gain obtained by combining the received signals at the two antenna branches depends on the amount of correlation between 
■ them, e.g., the lower the correlation the higher the diversity gain, Moreover, when the correlation coefficient is 1 or —1 it is 
Q\ ' said that the two signals are perfectly correlated, which implies that there is no diversity gain. In this case the joint probability 
• distribution of the signals is not well-defined and, to the best knowledge of the authors, a thorough discussion of these limit 
| cases is not found in the literature. Here, we fill this gap. 

We therefore investigate in this letter the representation of the joint probability distribution in the limit case of perfectly 
' correlated random variables. And show that this is achieved by Dirac's (5-function, [2 1. Based on this, we propose a representation 
of Dirac's <5-function as the ratio between bivariate distributions and the marginal distribution in the limit p — > ±1. An example 
of the evaluation of the limit is provided for the bivariate Rice distribution. 
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II. Dirac's (5-function representation 

Let X e R and Y £ R be two random variables defined on some probability space. Thus, the Pearson product-moment 
correlation coefficient between them is defined as, 101 

ctxvy 

where E[.] denotes the expected value operator, px = E[x] and p,y = E[x] denote the mean values, and ax and cry denote 
the standard deviations, i.e. 

ox = VE[x 2 } - nx 2 and a Y = ^E[y 2 }-p Y 2 . (2) 

Definition (Q]i ensures that — 1 < p < 1. 

Let fx.Y(x,y) denote the joint probability density functiorQof the two random variables X and Y. By definition we have 
that 

fx,r{x,y) = f Y \x(y I x)fx(x) = f x \y(x | y)fv(y), (3) 

where Jy\x{v I x) and fx\y(% I y) are the conditional distributions, fx(x) and fy(y) denote the marginal distributions of 
variables X and Y, respectively. 
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Let's further assume that random variables X and Y are identically distributed, with mean values of equal or opposite signs. 
We state then that the joint probability distribution fx,y(x, y) can be compactly represented as follows 

, - / fx(x)S(x T y) if P = ±1 and p x = ±p,y, (A , 
fx*( x >v)-\ if p = ±l and M x=TM^0, (4) 

where the upper and the lower signs should be evaluated separately, S(.) is the 5-function or Dirac's <5-function and fx(x)S(x=f 
y) = fr{y)5{yTx). 

Indeed, for perfectly correlated random variables p = ±1 we have from (Q]i that 

±a 2 = E[xy}TV 2 , (5) 

where the (+) and (— ) signs on the left hand side of correspond to p = 1 and p = — 1, respectively. On the righthand side 
the (+) and (— ) signs correspond to opposite sign mean values (p = px = —py) and equal sign mean values (p = px = Hy), 
respectively. 

From (01 we have that 

a 2 = E[x 2 ] -p 2 . (6) 

Consider now E[(x — ky) 2 ] = 0. Thus, solving for k and combining ([5]) and (|6]l gives 

E[xy] ± y/E 2 [xy]- E 2 [x 2 ] 
~ E\^~] 



(7a) 



±a 2 ±p 2 ± ^(±g ± p 2 ) 2 - (a 2 + p7f nM 
o A + p z 

where we have used our assumption of identically distributed variables. A straightforward analysis of (ITBl leads to three 
possibilities: (I), k = 1 if p = 1 and p = px = Py = or p = px = Py 7^ 0, (II), k = — lifp = — 1 and p = px = py = 
or p = px = —py 7^ 0, and (III), there is no (real) solution if p = 1 and p = px = —py 7^ or if p = —1 and 
p = px = Py 7^ 0. Hence, summarizing, we find that 

±1 if p = ±1 and p = px = ±Py, 

{0} if p = ±1 and p = p x = =fp Y ^ 0. ( ' 

Clearly, k = ±1 in the first case implies a linear, deterministic relationship between random variables X and Y with probability 
one. Thus, we can write 

fr\x(y I x) = S(y T x) and f x \ Y (x \ y) = S(x =p y), (9) 

since S(y ± x) = S(x ± y). Now, inserting (J8j into (0 gives 

fx,y(x,y) = f x (x)S(yTx) = f Y (y)$(xTy)- (10) 

On the other hand, k = {0} implies in the second case that there exist no linear, deterministic relationship between random 
variables X and Y and therefore a perfect correlation between the variables cannot be observed. Hence, in this case 

fx, Y (x,y) = 0. (11) 

This ends our derivation. 

Let us consider bivariate distributions of identically distributed random variables (with mean values of equal or opposite 
signs) as previously. We further assume they also depend on the correlation coefficient p as a parameter, and are well-defined 
functions for < \p\ < 1. However, such distributions are in general not well-defined for \p\ = 1 since terms of the form 
(l±p) _1 , (1 ± p)~ x / 2 , (1 ± p 2 )^ 1 and/or (1 ± p 2 ) -1 / 2 may appear in their arguments. Examples of such bivariate distributions 
are the Gaussian (for both complex and real variables) and the family of distributions derived from the complex Gaussian 
multivariate distribution , J3], such as the bivariate Rice , |4), the bivariate Nakagami, f5), and the bivariate Weibull distributions 
, 0. Henceforth, we denote such a bivariate distribution as Jx,y(x, y; p) to make a distinction between the joint distribution 
and bivariate distribution, i.e., we explicitly point out the dependence on p for < \p\ < 1. Therefore, we can write 

fxA^v) = fx, Y (x,y;p) if < \p\ < 1. (12) 

However, as we have shown above the singularities can be represented through the (5-function when p = 1 and px = Py, or 
p = — 1 and px = —py or p = ±1 and px = Py = 0. Thus 

fx.Y(x,y) = \im f x .Y(x,y;p) = fx{x)S(xTy)- (13) 
p— >±i 

Hence, Dirac's ^-function has the following representation 

xi -r- ^ v fx,Y{x,y;p) 

S(xTy)= lim — f / , n ■ (14) 
p^±! fx(±y) 
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Assume now that X c e C and Y c £ C are two random variables defined on some probability space. Thus, the Pearson 
product-moment correlation coefficient between them is defined as, (3) 

E[(x c -nx c ){yt -(J>y c )\ 

Pc = — , (15) 

&x c &Y c 

where E[] denotes the expected value operator, (.)* denotes complex-conjugate, px c — E[x c ) and py c = E[x c ] denote the 
mean values, and ax c and oy c denote the standard deviations, i.e. 

ox, = \/E[\x c - px c \ 2 } and a Yc = V E l\Vc - Py c \ 2 }- (16) 

We further assume that p c = p*, i.e., p c S R. In this case if X c and Y c are identically distributed with mean values of equal 
or opposite signs, then the following holds true 

\Pc\ = 1 p = ±1 if x = \x c \, y = ±\y c \, (17) 

where, the upper and the lower signs should be evaluated separately. The converse, i.e., \p\ = 1 \p c \ = 1, doesn't hold in 
general. This follows from the definition of the correlation coefficients ([1) and ( fTSV 

III. An example: the bivariate Rice distribution 
Consider the bivariate Rice distribution of two identically distributed random variables X <G R and 7el 

(1 a. K\ 2 th 2K (i+-K)(x 2 +» 2 ) 

fx,Y(x,y;p c ) = 9 l 2, e " 1+Pc " 2 " {1 - Pi> (18) 

2ir/3 2 {l - pi) 



2k 



J ^~ l ^ 2K(l + K){x>+ y 2 + 2xy cosg) j dg> 

2 2 

where /3 = %■ = ^- is the average power of, e.g., the signal received at one antenna, A' is the power of the line-of-sight(LOS) 
component to the scattered power and p c is the correlation coefficient J151 . and Iq(.) is the modified Bessel function of Oth 
order. We see from ([18} that fx,Y(x,y;p c ) is not well-defined for p c = ±1, so it should be evaluate in the limit p c — >• 1 to 
obtain a value if it exists. 

Evaluating the limit by regrouping terms and using the Algebraic Limit Theorem gives 

(1 4- K) TUP 2p 

lim f x , Y (x,y;p c ) = ( + ) l R2 G 1 (x,y) (19) 

2tt 



* j 0>(x ,,, w ^ £0±«2Cg^tjgs g fl j Jti 

where 

e 2/3(i- P ^) 

Gi(x,y) = lim -pr: n\i72' (20) 

p c (l+Jf)lj/(l-COB 8) 

e «RT 

G 2 (x,y,0) = lim — . , (21) 

The limits in d20l ) and d2Tb are of the form jj, which can be resolved by noticing that Dirac's-<5 can be represented in terms 
of the limit 

5(z) = lim^U. (22) 

Thus, 



Gi(a:, y) = 8{y/(l + K)/(2/3)(x - y)) , (23) 

and 



G 2 (x, y, 9) = S(y/xy(l + A)(l - cos B)/0) . (24) 
Let's first consider (|23l , which by applying the following property of Dirac's ^-function 

= E %jf' (25) 

n I da; I 
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where the x n satisfy g(x n ) — 0, is reduced to 

Gi(x,y) 

Similarly, applying d25T l to (l24l i we identify 



2(3 



:S(x - y), 



9(0) = y/((l + K)xy(l-cosO))/l3, 
with the modulus of the first derivative given by 

dg(6) 



(1(9 



(1 + K)xy 



23 



cos ■ 



(26) 



(27) 



(28) 



where we have used some standard mathematical identities and operations. Further, it is straightforward to see that the zeros 
of d27l i are given by 9 n = 2nn. Thus, 



G 2 (x,y,6) 

The shifting property of Dirac's ^-function 



2/3 



(1 + K)xy^ 



f(x)6(x - x )dx = f(x ), 



(29) 



(30) 



where T> is the domain of integration and xq € T>. Thus, inserting (l26l l and d29l into ( TT8l and applying (f30b we obtain the 
final result 



lim fx,v(x, y,p c ) = S(x - y) 

x(l + K) K (i+k)» 2 T / 
x — ^— — -e " 2,3 / x . 



(31) 



2if (1 + K) 



IV. Conclusions 



The Pearson product-moment correlation coefficient p plays a fundamental role in the analysis of the performance of diversity 
systems for wireless communications. We have shown that for bivariate distributions belonging to the Gaussian family, e.g., 
Rayleigh, Ricean, Nakagami and Weibull, and perfectly correlated random variables they can be represented by Dirac's 5- 
function if the corresponding linear dependence exists. The presented result can be used in the analytical and numerical analysis 
of the diversity performance of perfectly correlated diversity branches. In addition, for nearly perfectly correlated results can 
be obtained by expanding the joint probability distribution around p pa 1 or p m — 1. Finally, we believe that for the sake 
of completeness the limiting cases corresponding to p = 1 and/or p = — 1 should also be given when providing the joint 
probability distribution of two random variables. 
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